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Abstract 

\ The dispersive two-photon contribution to the Kl - ► decay am- 

plitude is analyzed, using chiral perturbation theory techniques and large- 
Nq considerations. A consistent description of the decays ir° — ► e + e~, 



fN| ■ n — * /J, + n and Kl — > /U + p is obtained. As a byproduct, one predicts 

>'. Br (7/ -» e+e") = (5.8±0.2) x 10 -9 and Br (if L -» e+e~) = (9.0±0.4) x 1(T 12 . 
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The rare decay Kl — > p + p~ has deserved a significant theoretical interest during the 
last three decades. It represents a potentially important channel to study the weak interac- 
tion within the Standard Model (SM), as well as possible effects of new physics, mainly in 
connection with flavour-changing neutral currents and CP violation. 

This decay proceeds through two distinct mechanisms: a long-distance contribution from 
the 27 intermediate state and a short-distance part, which in the SM arises from one-loop 
diagrams (W boxes, Z penguins) involving the weak gauge bosons. Since the short-distance 
amplitude is sensitive to the presence of a virtual top quark, it could be used to improve our 
present knowledge on the quark-mixing factor Vtd', moreover, it offers a window into new- 
physics phenomena. This possibility has renewed the interest in the study of the Kl — > p + p~ 
process in the last years. 

The short -distance SM amplitude is well-known |J. Including QCD corrections at the 
next-to-leading logarithm order ]2[] , it implies [|3| : 

Br (K L „♦„-)» = 0.9 x 10- (po - Pf (= j ^ j , (1) 

where po ~ 1-2 and p = p(l — A 2 /2), with p and A the usual quark-mixing parameters, 
in the Wolfenstein parameterization Q. The deviation of po from 1 is due to the charm 
contribution. Using the presently allowed ranges for mt and the quark-mixing factors, one 
gets [f| Br {Kl — > p + p~)sr> = (1.2 ± 0.6) x 10~ 9 . If this number is compared with the 
measured rate 0] 
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Br (K L -> = (7.2 ± 0.5) x 10~ 9 , (2) 

it is seen that the decay process is strongly dominated by the long-distance amplitude. 

Clearly, in order to extract useful information about the short-distance dynamics it is 
first necessary to have an accurate (and reliable) determination of the Kl —>■ 7*7* — > H + fi~ 
contribution. 

It is convenient to consider the normalized ratios 



where (3 = Jl — Amf /Mp. The on-shell 27 intermediate state generates the absorptive 
contribution 

ta[F(i>^r)] = i]n(£!). (4) 

Using the measured branching ratio 0, Br (K L — > 77) = (5.92 ± 0.15) x 10~ 4 , this implies 
the so-called unitarity bound: 

Br (K L -> n+fT) > Br (K L -> ^/Oa^ = (7.07 ± 0.18) x 10~ 9 . (5) 

Comparing this result with the experimental value in Eq. ([|), we see that Br (Kl — ► /x + /i~) 
is almost saturated by the absorptive contribution. 

One immediate question is whether the small room left for the dispersive contribution, 
Br (Kl — > /U + /i~)Dis = (0.1 ±0.5) x 10 -9 , can be understood dynamically. Naively, one would 
expect a larger value just from the intermediate 27 mechanism. This has motivated some 
recent speculations || about a possible cancellation between the long- and short-distance 
dispersive amplitudes, which could allow for additional new-physics contributions at short 
distances. 

The obvious theoretical framework to perform a well-defined analysis of the long-distance 
part is chiral perturbation theory (ChPT). Unfortunately, the chiral symmetry constraints 
are not powerful enough to make an accurate determination of the dispersive contribution 
PHill. The problem can be easily understood by looking at the Kl — > 77 amplitude, 

A(K L - 77) = c(ql q 2 2 ) e lfl e 2l/ q lp q 2a , (6) 

which, at lowest-order in momenta, proceeds through the chain Kl — > vr°, rj, rj' — > 27. The 
lowest-order — 0(p 4 ) — chiral prediction, can only generate a constant form factor c(qf, q%); 
it thus corresponds to the decay into on-shell photons (qf = q\ = 0) |fL2| : 

2G 8 afn 
7rd - rl 



c (°' °) = — c ^d , (7) 



c r ed = 1 - ^ 2 r *\ (c e - 2V2s e )(c e + 2V2 Pn s e ) + ( / 2 ^ (2V2c e + s e )(2V2p n c e - s e ) 
o{r v l) 6{r v , L) 
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where rp = Mp/M^ L , eg = cos8p and sg = sin#p, with dp pa —20° the rt-rj mixing angle. 
The global parameter G$ = 2~ 1 / 2 GFV ud V* s g$ characterizes JT3J the strength of the weak 
AS = 1 transition Kl — > tt°, 77, rj . 

In Eq. (0) we have factored out the contribution of the pion pole, which normalizes the 
dimensionless reduced amplitude c re d- The second and third terms in c re d correspond to the 
rj and rj contributions respectively. Nonet symmetry (which is exact in the large-iVc limit) 
has been assumed in the electromagnetic 27 vertices; this is known to provide a quite good 
description of the anomalous P — > 27 decays (P = 7r ,T],T]'). Possible deviations of nonet 
symmetry in the non-leptonic weak vertex are parameterized through p n 7^ 1. 

In the standard SU(3)l <%> SU(3)r ChPT, the rj contribution is absent and 6 P = 0; 
therefore, c ro d oc {3M 2 + M 2 — 4M|-), which vanishes owing to the Gell-Mann-Okubo mass 
relation. The physical Kl —>■ 77 amplitude is then a higher-order — 0(p 6 ) — effect in the 
chiral counting, which makes difficult to perform a reliable calculation. 

The situation is very different if one uses instead a U{3)l ® U(3)r effective theory p4| , 
including the singlet r\\ field. The large mass of the rj originates in the U(1)a anomaly 
which, although formally of 0(1/Nc), is numerically important. Thus, it makes sense to 
perform a combined chiral expansion [|l5j in powers of momenta and 1/Nc, around the nonet- 
symmetry limit, but keeping the anomaly contribution (i.e. the 7/1 mass) together with the 
lowest-order term. In fact, the usual successful description of the 77/7/ — > 2j decays fl6 



corresponds to the lowest-order contribution within this framework, plus some amount of 
symmetry breaking through 7^ f v > 7^ f n . The mixing between the rj 8 and rji states provides 
a large enhancement of the rj — > 27 amplitude, which is clearly needed to understand the 
data. In the standard SU (3)l<S>SU (3)r framework, the rj is integrated out and its effects are 
hidden in higher-order local couplings |17|]; since the 771 and r]s fields share the same isospin 
and charge, the singlet pseudoscalar does affect the r, dynamics in a significant way, which 
is reflected in the presence of important higher-order corrections ||18|| . These corrections are 
more efficiently taken into account within the U(3)l <S> U(3)r frameworkQ 

Taking sg = —1/3 {6 pa —19.5°), the ??-pole contribution in Eq. (0) is proportional to 
(1 — p n ) and vanishes in the nonet-symmetry limit; the large and positive rj contribution 
results then in c re d = 1.80 for p n = 1. With < p n < 1, the rj and rj contributions interfere 
destructively and c re d is dominated by the pion pole. One would get c rc d — 1 for p n ~ 3/4. 

The measured Kl — > 77 decay rate corresponds to |c(0,0)| = (3.51 ± 0.05) x 
10~ 9 GeV" 1 . With \G 8 \ = 9.1 x 10" 6 GeV" 2 , obtained from the 0(p 2 ) analysis of K -> 2vr 
p~3|1 , this implies cj^ = (0.84 ± 0.11). However, the fitted value of \G$\ gets reduced by 
about a 30% when 0{p A ) corrections to the K — > 2tt amplitudes are taken into account ||20| . 
This sizeable shift results mainly from the constructive tttt rescattering contribution, which 
is obviously absent in — > 77. Thus, we should rather use the corrected (smaller) \G$\ 



1 The contributions of the singlet pseudoscalar are particularly important in radiative transitions, 
owing to the presence of the rj exchange pole. A different situation occurs in the decays rj — > 
3tt, where the 0(p 4 ) corrections induced by 77—77' mixing are related to the pseudoscalar mass 
spectrum, and cancel to some extent with other contributions associated with the exchange of 
scalar particles [19|. 
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determination, which leads to c re ^ — (1.19 ±0.16). 

Leaving aside numerical details, we can safely conclude that the physical Kl — > 77 
amplitude, with on-shell photons, is indeed dominated by the pion pole (c re d ~ 1)- Although 
the exact numerical prediction is sensitive to several small corrections (p n 7^ 1, f w 7^ 
fv 7^ frfi s d ~ 1/3) and therefore is quite uncertain, the needed cancellation between the rj 
and 7]' contributions arises in a natural way and can be fitted easily with a reasonable choice 
of symmetry-breaking parameters. 

The description of the Kl — > 77 transition with off-shell photons is a priori more com- 
plicated because the q\ 2 dependence of the form factor originates from higher-order terms 
in the chiral lagrangian. This is the reason why only model-dependent estimates of the 
dispersive Kl — > transition amplitude have been obtained so far. At lowest-order in 
momenta, c(qf, gf) = c(0, 0) ; thus, the (divergent) photon loop can be explicitly calculated 
up to a global normalization, which is determined by the known absorptive contribution 
(i.e. by the experimental value of c(0, 0)). The model-dependence appears in the local con- 
tributions from direct KlI + 1~ terms in the chiral lagrangian P,|TTf (allowed by symmetry 
considerations), which reabsorb the loop divergence. 

It would be useful to have a reliable determination in some symmetry limit. The large- 
N c description of K L — > 7*7* provides such a possibility. At leading order, this process 
occurs through the 7r°, 77,77' poles, as represented in Fig. 1. Therefore, the problematic 
electromagnetic loop in Fig. 1(a) is actually the same governing the decays tt° — > e + e~ and 
77 — > and the unknown local contribution (Fig. 1(b)) can be fixed from the measured 



rates for these transitions p2|,p3|. In fact, the same combination of local chiral couplings 
shows up in both decays |22[, leading to a relation that is well satisfied by the data. 

Although the 7/ — > transition introduces additional chiral couplings, they are sup- 
pressed by at least one more power of 1/Nq- Thus, in the large-iV^ limit the different 



electromagnetic P — > 11 decays get related through the same counterterms [22 



£c.t. = |^ (h^l) {xi Tr (Q 2 {U\ 8,U}) + X 2 Tr (QU^Qd.U - Qd^Qll) } , (8) 

where Q = diag(2/3, — 1/3, — 1/3) is the quark electromagnetic charge matrix and U = 
exp ^i\/2$//^ the usual U(3)l® U(3)r matrix describing the pseudoscalar nonet. 

Nonet symmetry should provide a good estimate of the ratio R(Kl — > Since 
Kl — > 77 is dominated by the pion pole, we can reasonably expect that symmetry-breaking 
corrections would play a rather small role. In any case, this symmetry limit allows us to 
investigate whether the tiny dispersive contribution allowed by the data is what should be 
expected from the 27 intermediate state. 

In this limit, all R(P — > ratios are governed by the same dispersive amplitude^: 



Re[F(P^l + r)\ = -=-ln 2 + j-Li 2 ^— ^ + -7— + 31n — + x(jj) , (9) 





2 Notice that our result differs slightly from those quoted in Refs. |23| and |24j . We agree with the 
numerical expression given in Ref. [p2]. 
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where xil 1 ) = — (XiO") + X2M + 14) /4 is the relevant local contribution, with Xiit 1 ) (i=l,2) 
the corresponding chiral couplings renormalized in the MS scheme. The /1 dependence of 
the x{p) an d In (m;//i) terms compensate each other, so that the total amplitude is fi- 
independent. 

Table | shows the fitted values of x(M p ) from the three measured ratios R(tc° — > e + e~), 
R(r] — > and R(Kl — > /i + /i~). Subtracting the known absorptive contribution, the 

experimental data provide two possible solutions for each ratio; they correspond to a total 
positive (solution 1) or negative (solution 2) dispersive amplitude. We see from the Table 
that the second solution from the decay n° — > e + e~ is clearly ruled out; owing to the 
smallness of the electron mass, the logarithmic loop contribution dominates the dispersive 
amplitude, which has then a definite positive sign (an unnaturally large and negative value 
of x{M p ) is needed to make it negative). The large experimental errors do not allow to 
discard at this point any of the other solutions: the remaining value from tt° — > e + e~ is 
consistent with the results from the rj — > ji + /i~ and K L — > fi + fi~ decays, and these are also 
in agreement with each other if the same solution (either the first or the second) is taken for 
both. We see that, in any case, the three experimental ratios are well described by a common 
value of x{M p ). In this way, the experimentally observed small dispersive contribution to the 
K L — > decay rate fits perfectly well within the large-iVc description of this process. 

We have not considered up to now the short-distance contribution to the K L — > n + fi~ 
decay amplitude ||. This can be done through a shift of the effective x(M p ) value|: 

X{M P U = x(M p ) - 5xsd , (10) 
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^Xsd w 1.7 (p - p) 

For the allowed range \p\ < 0.3, one has Sxsd ~ 1.8 ± 0.6, which allows to exclude the 
solution 2 for x(M p ) obtained from 77 — > fi + ^. The solution 1, on the contrary, is found to 
be compatible with the results from K L — > p + /i~, and can be used to get a constraint for 
&XSD- Indeed, taking as the best determination 

X (M p ) = 5.5^, (11) 

the first solution for K L — ► leads to 

<JXSD = 2.2lJi, (12) 

in agreement with the 5xsd value quoted above. The second solution for K L — > p + p>~ 
appears to be less favoured, yielding <5xsd = 3.6 ± 1.2; this shows a discrepancy of about 
1.4(7 with the short-distance estimate. Notice that the precision of the result in (|I"2"D is still 
relatively low. However, the errors could be reduced by improving the measurements of the 
f] — > and Kl — > fi + fi~ branching ratios. 



5 The relative sign between the short- and long-distance dispersive amplitudes is fixed by the 



known positive sign of <?§ in the large-iV^ limit [ 25 1 . 
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Finally, once the local contribution to the P — > l + l decay amplitude has been fixed, it 
is possible to obtain definite predictions for the decays into e + e~ pairs: 

Br(7r° -> e+e") = (8.3 ± 0.4) x lO" 8 , Br(^ -> e+e") = (5.8 ± 0.2) x lO" 9 , , , 
Br(K L -> e+e") = (9.0 ± 0.4) x lO" 12 . [ 6) 

In the same way, the amplitudes corresponding to the rj' decays are found to be Br(r/ — > 
e+ e -) = (1.5±0.1)xlO- 10 andBr(?/ = (2.1 ±0.3) x lO" 7 . However, in view of the 

large mass of the rj', these predictions could receive important corrections from higher-order 
terms in the chiral lagrangian. 

To summarize, in the nonet symmetry limit it is possible to make a reliable determination 
of the ratios R(P — > l + l~), at lowest non-trivial order in the chiral expansion. A consistent 
picture of all measured P — > l + l~ decay modes is obtained, within the SM. The present data 
allow to pin down the size of the relevant chiral counterterm and to get a constraint on the 
short-distance contribution to the Kl — > amplitude. However, this constraint is found 

to be rather weak; more precise measurements of the 77 — > fi + fi~ and Kl — > fi + fi~ branching 
ratios would be necessary in order to improve the bounds for 5xsd obtained through the SM 
box and penguin computations. In addition, a more detailed investigation of the theoretical 
uncertainties is needed to quantify how precisely the short-distance Kl — * fi + fi~ amplitude 
can be inferred from the data. Although it seems difficult to achieve a theoretical precision 
good enough for making useful tests of the SM flavour-mixing structure, it is worth to try. 
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FIGURES 




(b) 

FIG. 1. (a) Photon loop and (b) associated counterterm contributions to the Kl 
process. 
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TABLES 



TABLE I. Fitted values of xiMp) from different R(P — > l + l ) ratios. The numbers quoted for 
Kl — ► refer to the difference %(M p ) — #xsd- 





x(M p ) [Solution 1] 


x(M p ) [Solution 2] 


vr° - e+e" 


41* 

c c +0.8 
°-° -1.0 
q q+0.9 
* -0.7 


-24 ± 5 
-0 8 +L0 

u -° -0.8 
1 q+0.7 
1V -0.9 
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